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Abstract 

(N ' 

J> ' We present, using the statistical model, an analysis of the production of light nuclei, 

hypernuclei and their antiparticles in central collisions of heavy nuclei. Based on 
. these studies we provide predictions for the production yields of multiply-strange 

CN \ light nuclei. 

d ' 



1 Introduction 
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One of the major goals of ultrarelativistic nuclear collision studies is to obtain informa- 
tion on the QCD phase diagram [TJ. Currently, one of the most direct approaches is the 
investigation of hadron production. Hadron yields measured in central heavy ion collisions 
from AGS up to RHIC energies can be described very well [2|3Pf5|7f8|9|10|ll|IT2] within 
a hadro-chemical equilibrium model. In our approach p|¥|7ffT2|13|14] the only parameters 
are the chemical freeze-out temperature T and the baryo-chemical potential //& (and the 
fireball volume V, in case yields rather than ratios of yields are fitted). Other approaches 
[51l^|ll|15|ITrj] employ (several) other, non-thermal, parameters. For a review see [T7] . 

The main result of these investigations was that the extracted temperature values rise 
rather sharply from low energies on towards ^s^n —10 GeV and reach afterwards con- 
stant values near T=160 MeV, while the baryochemical potential decreases smoothly as 
a function of energy. This limiting temperature [18] behavior suggests a connection to the 
phase boundary and it was, indeed, argued [TH] that the quark-hadron phase transition 
drives the equilibration dynamically, at least for SPS energies and above. For the lower 
energies, the quarkyonic state of matter [20] could complement this picture by providing 
a new phase boundary at large fib values. The conjecture of the tricritical point [21] was 
put forward in this context. 



The importance of measurements at very high energies to obtain information on the exis- 
tence of a limiting temperature of excited hadronic matter produced in nuclear collisions 
was pointed out early [22|23f24|l2o'] based on analysis of particle spectra at the Bevalac 
(see also the review [26]), from pions to heavier complex nuclei. 

At first glance, it may seem inappropriate to use the chemical freeze-out concept for 
light nuclei, as their binding energies are a few MeV, much less than the chemical freeze- 
out temperatures of 100-170 MeV. We note, however, that the relative yield of particles 
composed of nucleons is determined by the entropy per baryon, which is fixed at chemical 
freeze-out. This has been first recognized already 30 years back [22] and was subsequently 
further substantiated in [23], constituting the basis of thermal analyses of yields of light 
nuclei [2~Tf28] . It is entropy conservation, and not the difference between the binding energy 
and temperature of the system, which governs the production yields in this case After 
chemical freeze-out, entropy is conserved. 

It was also noted then that the yields obtained within the thermal model are in close 
agreement to those from coalescence models [271129"] . The thermal model studies were 
already at that time extended to nuclei carrying strangeness (hyperons in replacement 
of nucleons) and even hypothetical objects with roughly equal number of up, down and 
strange quarks (strangelets). At the same time, a vigorous line of theoretical investigations 
on the existence of multi-strange hypernuclei, or MEMOs [30|31|32|l33j was established. 

Recently, the first measurement of the lightest (anti)hypernucleus, (anti) hyper-tritium, in 
high-energy nucleus-nucleus collisions was achieved by the STAR experiment at the RHIC 
|34j . This measurement opens up a very interesting new regime for tests of particle produc- 
tion at chemical equilibrium. At relatively low beam energies, where the baryo-chemical 
potential and, hence, the baryon density is maximum (FAIR energy regime) objects with 
a large number of baryons and moderate strangeness may be abundantly produced [33] . 
At RHIC and LHC energies production of objects with moderate (anti)baryon number 
and large strangeness content may be expected. In this paper we investigate and predict 
within the thermal model the production yields of heavy baryons and anti-baryons and 
in particular of hypernuclei and their antiparticles and confront these calculations with 
all presently available data ranging from AGS to RHIC energies. 



2 Preliminaries 

The measurement of the production yields of light nuclei (and anti-nuclei) without 
strangeness in central nuclear collisions provides significant constraints on thermal model 
parameters, in particular on the value of the baryo-chemical potential /if,. This is most eas- 
ily seen when one recognizes that yield ratios such as n He/ n He scale like exp [— (2n/x&/T)]. 
In Fig. [1] we show the updated thermal fit to the hadron yield data measured at RHIC 
(v /s A r A r= 200 GeV) including the newly-measured [34] yield ratio 3 He/ 3 He. Including this 
ratio significantly narrows the range of possible fib values, while T and V of the new fit 
remain unchanged (T=164 MeV, V=1960 fm 3 ) compared to our earlier fit |13j . Quanti- 
tatively, the new fit leads to /i;>=24±2 MeV, while without the ratio 3 He/ 3 He, /if,=30±4 
MeV [13] . The quality of the present fit is similar to that of the earlier one (which had 
X 2 /dof=29.7/12). This result supports previous findings at lower energies [27|28|29j . We 
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Figure 1. Hadron yields in comparison with the thermal model fit of combined data (excluding 
K*, £*, A*), for the RHIC energy of ^/snn =200 GeV. The ratio 3 He/ 3 He, recently measured 
by the STAR experiment [M] is included in the fit. 



stress that the agreement between the experimental value and the calculated one for the 
ratio 3 He/ 3 He is a powerful argument that indeed entropy conservation governs the pro- 
duction of light nuclei. If one were to use a temperature comparable to the binding energy 
per nucleon, that is T=5 MeV, the calculated ratio would be 3.1-10" 13 , while it is 0.415 
for T=164 MeV, see Fig. [TJ 
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Figure 2. Energy dependence of anti-baryon to baryon yield ratios. The lines are thermal model 
results as described in the text. The symbols represent measured data. 

In Fig. [2] we show that predictions using the thermal model can be used to describe 
quantitatively the measured energy dependence oip/p, d/d, and A/A yield ratios over a 



very wide energy range The calculations, here and in the following, are performed using 
the parametrizations for T and jib established in [13] based on fits of midrapidity data 
in central collisions. The penalty against anti-particles at lower energies, well described 
within the thermal approach, is also drastically exhibited in this figure. Note, in particular, 
the very good agreement between the model and the measurements at AGS [35] for the d/d 
ratio, which extends the the range of the model agreement over 8 orders or magnitude (for 
the other measurements see references in [T2][T3"] ). It is thus natural to extend such analyses 
to light nuclei containing strangeness. Note that, at the lowest energies, the canonical 
strangeness suppression is important and is incorporated in our model as described in 
[12]. In the following we will therefore use the thermal model with the parameters as 
discussed to analyze the production of hypernuclei and their anti-particles and confront 
model predictions with the now available data. 



3 Production of Hypernuclei at RHIC Energy 

In Table [1] we show a comparison of the measured data [31] and model calculations 
for the best fit thermal parameters discussed above. The yield ratio \H /\H is as well 
reproduced by the model as the ratio 3 He/ 3 He. On the other hand, the measured ratios 
of (anti)hyper-tritium to (anti) 3 ife are larger than predicted in the model by about two 
standard deviations, using the statistical and systematic uncertainties quoted for the data. 
Table 1 

Ratios at RHIC energy, y / s]v r /v=200 GeV. The experimental values are from the STAR exper- 
iment |34j and contain statistical ans systematic errors. The errors for the model calculations 
correspond to the errors of the fit for the baryochemical potential, fib = 24 ± 2 MeV. 



Ratio 


Experiment 


Model 


3 He/ 3 He 


0.45±0.02±0.04 


0.42±0.03 


\H/\H 


0.49±0.18±0.07 


0.45±0.03 


\H/ 3 He 


0.82±0.16±0.12 


0.35±0.003 


\H/ 3 He 


0.89±0.28±0.13 


0.37±0.003 



To shed more light on the situation we turn now to the energy dependence of (strange) 
baryon production. In Fig. [3] we show the experimental energy dependence of the A/p 
and d/p ratios and confront these data with our thermal model predictions. The degree 
of agreement between data and calculations is impressive. We also include in this figure 
thermal model predictions for the energy dependence of the ratio \H/ 3 He and \H j 3 He. 
The broad maximum around ^/snn —5 GeV for the ratio \H / 3 He has the same origin 
as the maximum in the K + /ir + ratio, namely it arises as a consequence of strangeness 
neutrality condition, imposed in our model, and a competition between rising T and 
decreasing fib [13]. We also note that a slightly less prominent maximum is likely to 
survive even if one relaxes the condition of strangeness neutrality, as demonstrated in 
|33j . At high energies the value for the two ratios approach each other, as expected for 
decreasing values of fib at a nearly constant temperature. 

In Fig. H] we show the measured energy dependence for the 3 He/ 3 H and the 
\H I '( 3 H e(A / 'p)) ratio. This double ratio was suggested by the authors of [M] in the ex- 
pectation that dividing out the strange to non-strange baryon yield should result in a 
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Figure 3. Energy dependence of various baryon yield ratios. The lines are calculations, the 
symbols are experimental data. 
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Figure 4. Energy dependence of nuclei and hypernuclei production ratios. The data points are 
extracted from ref. [31], the lines are our model calculations. Note the discrepancy between data 
and model for the double ratio involving hyper-tritium, where the continuous line represents the 
physical case, while the dashed line represents the case without the important contribution from 
feed-down from strong decays on the A/p ratio (see text). 

value near unity. The data are compared to thermal model predictions. Note that there is 
negligible feed-down from heavier states into states with baryon number 3. As expected, 
the measured energy dependence of the 3 He/ 3 H is well reproduced by the model calcu- 
lations. On the other hand, the discrepancy between thermal model predictions and data 



for the \H / ( 3 H e(A / p)) ratio is apparent (red line). It is important to realize that the 
ratio A/p is significantly influenced by feed-down from strong decays of excited baryonic 
states, leading to a value for the double ratio significantly below unity. The red line in 
Fig. H] contains such feed-down, the blue dashed line represents a calculation where the 
feeding is artificially left out. The feed-down from strong decays increases the A/p ratio 
and, hence, reduces the overall ratio. The STAR collaboration actually measures a ratio 
close to 1, above the thermal model prediction by twice the error quoted by the experi- 
ment. Interestingly, results from the E864 collaboration [36] (as shown in ref. [M]) at AGS 
energy are, albeit with large uncertainties, consistent with the thermal model prediction. 
The discrepancy at RHIC energy, if experimentally established, would point to a new 
production mechanism not contained in the thermal approach and not present at lower 
beam energies. The possible existence of an excited J 71 " = 3/2 + state of (anti)hyper-tritium 
has been recently pointed out to us [37]. This excited state could contribute via decay to 
the ground state, and would lead to close agreement between model and data. Further 
measurements at RHIC and, very soon, LHC energy are eagerly awaited to shed light on 
the situation. 




Figure 5. Energy dependence of hypernuclei to A yield ratios. 

To complete our studies we show, in Fig. [U predictions within the thermal model for the 
energy dependence of the production yield of multistrange light hypernuclei [27] relative 
to A hyperons. These ratios exhibit a pronounced maximum in the FAIR energy regime, 
which is the consequence of a competition between a strong increase (followed by satura- 
tion) of T and a strongly decreasing \i h (see also the discussion above). In addition, the 
canonical suppression, arising from the condition of local strangeness conservation, leads 
to reduced yields at low energies. In case of hyper-tritium production, there is no maxi- 
mum, since it is mainly determined by the strong energy dependence of fib at low energies. 
It is larger at the (low) FAIR energies by two to three orders of magnitude compared to 
RHIC and LHC energies. Even larger are the differences between low and high energies 
for the production of the exotic multi-hyperon states. 



In Fig. [6] we show, as a function of energy, predictions of yields at midrapidity per one 
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Figure 6. Energy dependence of predicted hypernuclei yields at midrapidity for 10 6 central 
collisions. The predicted yields of s He and ^He nuclei are included for comparison, along with 
the corresponding anti- nuclei (dashed lines). 

million central collisions. The volume at chemical freeze-out is that from our fits of yields 
|14j . At FAIR energy the production yields of exotic nuclei is maximal, although the 
absolute yields are still rather small. As an example, for \\^He, the rate of production 
at 10 6 central Pb+Pb collisions per second is about 60 per month, for a reasonable duty 
factor of the accelerator. Assuming a reconstruction efficiency of the order of a percent, 
this implies a few candidates per year of data taking, clearly at the edge of achievability. 

At the LHC, (anti-) 4 He and their corresponding hypernuclei are experimentally accessible. 
For the LHC energy of 2.76 TeV of the present data taking, we predict ratios 3 He/ 4 He and 
3 He/ 4 He of 2.76-10" 3 and 2.70T0 -3 , respectively (to be compared to the corresponding 
values for the RHIC energy of 200 GeV of 3.13T0 -3 and 2.37-10" 3 ). These predictions can 
also be used as guideline for expectation in pp collisions at LHC energy, where one could 
estimate, in the grand-canonical limit, yields reduced by a factor of the order of 200-400 
compared to Pb+Pb collisions. This is compensated by the much larger number of pp 
collisions (about 10 9 events) which can be inspected at LHC (for a running time of 10 7 
s per year of operation). This should allow, at the LHC, a measurement of the yields of 
produced (anti-)hypernuclei up to mass number 4 in pp and Pb-Pb collisions and provide 
a detailed test of our predictions. 




4 Conclusions 



We have demonstrated that the yield of light nuclei and their anti-particles are well 
reproduced with thermal model calculations employing parameters established from the 



analysis of general hadron production in relativistic nuclear collisions. As shown above, 
such ratios can be used to provide a precision constraint of the baryo-chemical potential 
fib- We have furthermore shown that the newly measured yield ratio \H /\H is also well 
described with the thermal approach, while the ratio \H/ 3 He which is reproduced at AGS 
energy is significantly underpredicted at RHIC energy. The origin of this discrepancy is 
currently not clear and needs further study. 

Our studies have also indicated interesting energy dependence in such yields and ratios. 
In particular, particles with large baryon number and moderate strangeness are produced 
in significant numbers at FAIR energy. 

The hyper-nuclei program, started by the STAR experiment at RHIC, has made these 
studies very topical. Although significant questions remain, it is clear that the study of the 
production of complex nuclei with and without strangeness in relativistic nuclear collisions 
can open a new chapter in the quest to understand the relation of particle production to 
the QCD phase boundary. The thermal model predictions can hopefully soon be tested 
also at the LHC energy with the data already collected in 2010. 
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